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Abstract. We obtain a complete asymptotic expansion of the integrated density of 
states of operators of the form H = (— A)™ + _B in M.'^. Here w > 0, and B belongs to a 
wide class of almost-periodic self-adjoint pseudo-differential operators of order less than 
2w. In particular, we obtain such an expansion for magnetic Schrodinger operators with 
either smooth periodic or generic almost-periodic coefficients. 



1. Introduction 

In [5], two of the authors of this paper have obtained the complete power asymptotic 
expansion of the integrated density of states of Schrodinger operators 

H = -A + V 

acting in M'^ assuming that the real-valued potential V is either smooth periodic, or 
generic quasi-periodic, or belongs to a reasonably wide class of almost-periodic functions 
(see [5] for a complete set of conditions on V as well as the previous history of the 
subject). The main aim of the current paper is to extend the results of [5] to a more 
general class of operators. We give a detailed description of this new class in the next 
section; here, we list the main properties of the operators belonging to it. 

(i) We consider perturbations of the Laplacian, or any positive power of the Laplacian. 
More precisely, we work with operators of the form 

where B is a differential or pseudo-differential operator of order k, < 2w. Here H is 
self-adjoint and belongs to the standard algebra of almost-periodic pseudo-differential 
operators, see e.g. [7] and [8]. 

(ii) If i? is a differential operator, we assume that its coefficients satisfy the same 
conditions the potential V had to satisfy in [5] (for example, the coefficients can be 
smooth periodic, or generic quasi-periodic functions). In particular, periodic magnetic 
Schrodinger operators are covered by our results. 
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(iii) If B is pseudo-differential, we assume tliat it is a classical pseudo-differential 
operator, or, more generally, the operator of classical type. By the latter we mean that 
the symbol of B admits an asymptotic decomposition in powers of |^| when |^| — )■ 00; 
however, these powers do not have to be integer. Note that operators with the relativistic 
kinetic energy ^^/(^^W~^-^K)^~+^rn^ are admissible for (almost-)periodic smooth A and 
m ^ 0. 

Under these assumptions we prove that the integrated density of states A^(A) has the 
complete asymptotic expansion (2.17). This expansion contains powers of A and powers 
of In A; the values of the exponents in the powers of A depend on the form of B, whereas 
logarithms are raised to integer powers smaller than d. Sometimes (as in the case of the 
magnetic Schrodinger operator) we can guarantee that the logarithmic terms are absent 
(i.e., the corresponding coefficients are zero). 

Remark 1.1. The main reason why we need assumption (iii) is to match asymptotic 
expansions in different intervals /„ in Section 3. If we did not have assumption (iii), we 
would have obtained the asymptotic expansions containing the general 'phase volumes' 
(hke in [9]), and it is not clear how to relate the expansions obtained in different intervals 

In. 

One immediate and shghtly unexpected corollary of (2.17) is as follows: 

Corollary 1.2. Suppose, H = (— A)^ -|- B with B being periodic and either differential, 
or pseudo- differential operator of classical type. Then for sufficiently large A the spectrum 
of H is purely absolutely continuous. 

Proof. Since H is periodic, the general Floquet-Bloch theory implies that the spectrum of 
H is absolutely continuous with the possible exception of eigenvalues of finite multiplicity. 
If A is such an eigenvalue, the integrated density of states has a jump at least |r^| at A, 
where is the lattice dual to the lattice of periods of H. Due to (2.17), this cannot 
happen for large A. □ 

The approach of our paper is similar to the one of [5] . In particular, we use the method 

of gauge transform developed in [9], [10], and [6]. Nevertheless, there are plenty of new 
(mostly technical, but sometimes ideological) difficulties arising because the operator B 
is no longer bounded and no longer local. One example of the new methods employed 
in this paper is the proof of Lemma 10.5: not only this proof works for unbounded S, it 
also makes Condition D from [5] redundant. The biggest increase in technical difficulties 
comes in Section 10 where we express the contribution to the density of states from 
various regions in the momentum space as certain complicated integrals and then try to 
compute these integrals. As a result, our paper is technically more complicated than [5] 
(which already was quite difficult to read). Thus, we have reluctantly abandoned the 
idea of making our paper completely self-contained; we will skip all parts of the argument 
which are identical (or close) to corresponding parts of [5] and refer the reader to that 
paper. Nevertheless, we will present all the definitions and properties of the important 
objects. 
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Remark 1.3. Throughout the article we employ the convention that, if some statement 
is given without a proof, then an analogous statement can be found in [5] , and the proof 
is the same up to obvious modifications. It comes withoTit saying that the reader is 
strongly encouraged to read the article [5] first, before attempting to read this paper. 
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2. Preliminaries 
For w > we consider the operator 

(2.1) H^{-A)'" + B 

acting in L2(M^). The action of the pseudo-differential operator B on functions from the 
Schwarz class S(R'') is defined by the formula 

(i^/)(x) := {27r)-'/'J 6(x,Oe^^-(3^/)(|)d|. 
Here 3^ is the Fourier transform 

(3^/)(|) := (27r)-'^/2 J e-$x^(x)rfx, ^ e 

the integration is over W^, and b is the symbol of B. We assume that 6(x, x, ^ e M.^, 
is a smooth almost-periodic in x complex-valued function and, moreover, that for some 
countable set © of frequencies we have 

(2.2) 6(x,|) = ^6(0,Oe,(x) 

where 

(2.3) e,(x) := e^«^ 
and 

b{0,^) :=M.(6(x,0e-«(x)) 

are the Fourier coefficients of b (here Mx is the mean of an almost-periodic function of x). 
We assume that the series (2.2) converges absolutely, and that b satisfies the symmetry 
condition 

b{0,^)^b{-d,^ + e), 
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SO that the operator B is formally sclf-adjoint. For i? > let 1^^ be the indicator 
function of the ball := : |^| < i?}. We assume that there exists a constant Cq 
such that 

and that 

(2.4) (1 - = E i^r^^(^'^/i^i)' 

where J is a discrete subset of (— oo, x] with 

(2.5) O^x <2w 

(the first inequality here is assumed for convenience without loss of generality), and 
6t(x, ry) are smooth functions on R*^ x S'^"^ almost-periodic with respect to x. 
Let 

(2.6) w :^ {w + x)/2. 
We introduce x e C°°(R+) so that 



(2.7) x{r) 



r, r ^ Co, 
0, r^Co/2. 



Remark 2.1. Increasing Co if necessary, we can guarantee that for any J C J and any 
C the operator B with the symbol b given by 

(2.8) &(x,0:=E(^(l^l))'E^^(^'^/l^l)^«(^) 
satisfies 

(2.9) (-A)* - \B\ ^ 0. 
We also assume that the coefficients in the expansion 

(2.10) b,{^,r]) = J2U^,v)M^), xeR^ T7eS'^-\ ceJ 

ee& 

can be represented by a series 

(2.11) 6,(0,771,...,%)= E^^HW---^? 

which converges absolutely in a ball of radius greater than one of M^. 

Under the above assumptions if is a selfadjoint operator on the Sobolev space H^"'(M'^). 
We are interested in the asymptotic behaviour of its integrated density of states N{X) 
as the spectral parameter A tends to infinity. 
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Definition 2.2. Let e(A;x, y) be the kernel of the spectral projection of H. We define 
the integrated density of states as 

N{\) :=M,(e(A;x,x)). 

It was proved in Theorem 4.1 of [8] that for differential operators this definition agrees 
with the traditional one (at least at its continuity points). The following lemma is proved 
at the end of Section 4 of [5]. 

Lemma 2.3. (a) IfA^B, then N{X; A) ^ N{\; B) . 

(b) Suppose, A = a(x, D) and U = m(x, D) are two pseudo-dijferential operators with 
almost-periodic coefficients. Let operator A be elliptic self- adjoint and operator 
U be unitary. Then N{X;A) = N{\] W^AU). 

Without loss of generality we assume that © (recall (2.2)) spans M"*, contains and 
is symmetric about 0; we also put 

(2.12) 0fe := + + --- + e 

(algebraic sum taken k times) and ©oo := ^k®k = ^(0)) where for a set 5 C by 
Z{S) we denote the set of all finite linear combinations of elements in S with integer 
coefficients. The set ©oo is countable and non-discrete (unless B is periodic). We will 
need 

Condition A. Suppose that di, . . . ,dii e ©oo- Then Z{di, . . . , 6^) is discrete. 

It is easy to see that this condition can be reformulated like this: suppose, 0i, . . . ,0d & 
©oo- Then either {Oj} are linearly independent, or Yl'j=i^j^j ~ 0' where rij G Z and 
not all Uj are zeros. This reformulation shows that Condition A is generic: indeed, if 
we are choosing frequencies of b one after the other, then on each step we have to avoid 
choosing a new frequency from a countable set of hyperplanes, and this is obviously a 
generic restriction. Condition A is obviously satisfied for periodic B, but it becomes 
meaningful if B is quasi-periodic (i.e., if it is a linear combination of finitely many 
exponentials) . 

If © and J are finite. Condition A is all we need. If, however, any (or both) of these 
sets is infinite, we need other conditions which describe, how well B can be approxi- 
mated by operators with quasi-periodic symbols. In the proof we are going to work with 
quasi-periodic approximations of B, and we need these conditions to make sure that all 
estimates in the proof are uniform with respect to these approximations. 

We introduce 

b,(6>) := sup |6,((9,r7)|, 6> e ©. 
|„|=i 



6 



S. MOROZOV, L. PARNOVSKI & R. SHTERENBERG 



Condition B. Let k be a positive integer. Then there exists Rq ^ Cq such that for each 
p > Rq there exist a finite symmetric set C (0 fl !B(p^/^)) (where S(r) is the ball of 
radius r centered at 0) and a finite subset J G J with 

(2.13) card J ^p^/'^ 

such that 

(2.14) Yl (1 + mT^'iRor^'Kie) < p~'. 

(e,t)e(0xj)\(0xj) 

The last condition we need is a version of the Diophantine condition on the frequencies 
of b. First, we need some definitions. We fix a natural number k (the choice of k will be 
determined later by the order of the remainder in the asymptotic expansion) and denote 

0j^ := 0^ \ {0} (see (2.12) for the notation). We say that QJ is a quasi-lattice subspace 
of dimension m, if QJ is a linear span of m linearly independent vectors 0i, . . . , 6^ from 
0^. Obviously, the zero space (which we will denote by X) is a quasi-lattice subspace of 
dimension 0, and is a quasi-Iatticc subspace of dimension d. We denote by V^, the 
collection of all quasi-lattice subspaccs of dimension m and put V := U^V^. If ^ G M'^ 
and QJ is a linear subspace of R*^, we denote by the orthogonal projection of ^ onto QJ, 
and put 5}-*- to be an orthogonal complement of QJ, so that ^ = +^2j-l- Let 2J, it G V. 
We say that these subspaces are strongly distinct, if neither of them is a subspace of 
the other one. This condition is equivalent to stating that if wc put Q2J := QJ fl il, then 
dimOU is strictly less than dimensions of QJ and 11. We put (p = 0(5J,il) G [0,7r/2] to 
be the angle between them, i.e. the angle between QJ OH and il W, where QJ 21J 
is the orthogonal complement of 22? in QJ. This angle is non-zero iff and Q2J are 
strongly distinct. We put s = s{p) = s(0^) := inf sin (0(QJ, it)), where the infimum 
is over all strongly distinct pairs of subspaces from V, R = R{p) := sup^^Q, \d\, and 
r = r{p) := inf^^~^ \d\. Obviously, 

(2.15) i?(p)=0(pV^), 
where the implied constant can depend on k and k. 

Condition C. For each fixed k and k the sets 0j^ can be chosen in such a way that for 
sufficiently large p the number of elements in 0j^ satisfies card 0j^ ^ p^/'^ and we have 

(2.16) s{p) ^ p-^l^ 
and 

r{p) > p-"\ 

where the implied constant (i.e. how large should p be) can depend on k and k. 

Remark 2.4. Note that Condition C is automatically satisfied for quasi-periodic and 
smooth periodic B; see [5] for further discussion of this condition. 
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Condition A implies the following statement, which will be used crucially in our con- 
structions. 

CoroUeiry 2.5. Suppose, 6i, . . . ,9i G 0^, / ^ d — 1. Let QJ be the span of di, . . . ,6i. 
Then each element of the set is a linear combination of 0i, . . . ,0i with rational 

coefficients. Since the set ©j^flQJ is finite, this implies that the set Z(0j^nQJ) is discrete 
and is, therefore, a lattice in QJ. 

From now on, we always assume that B satisfies all the conditions from this section; 
we will also denote 

p := X'/'^. 
Now we can formulate our main theorem. 

Theorem 2.6. Let H he an operator (2.1) satisfying Conditions A, B and C. Then for 
each K & there exists a finite positive integer L and a finite subset Jq d J such that 

(2.17) 

A^(p2"') 

d-l L [K+d+{2-2w)h+Li+-+th] 
q=0 h=0 i,i,...,i,h&Jo j=0 

as p ^ CO. 

Remark 2.7. The powers of p present in (2.17) are equal to d+{2—2w)h+Li + - ■ ■ + Lh—j, 
and the first impression is that there are far too many of them (indeed, a priori the set 
of all such powers can be dense in R, for instance). However, many of these powers are, 
in fact, spurious (i.e. the corresponding coefficients C^q^^j^ are zero). This happens, for 
example, when + (2 — 2w)h + t-i + ■ ■ ■ + Lh — j > d (for obvious reasons). Equally 
obviously, these powers do not 'multiply' when we increase K. This means that if 
Ki < K2, then expansion (2.17) with K = K2 does not contain extra terms with 
d + (2 — 2w)h + ii H \- I'h— 3 > —Ki, compared to this expansion for K — Ki. 

In the case of magnetic Schrodinger operators. Theorem 2.6 and calculations similar 
to those of [1] and [5] imply that most of the terms in (2.17) will indeed disappear: 

CoroUeiry 2.8. For each K eN we have: 

(2.18) A^(A) = A'^/^ (Ca + ^ ejX'^ + o(A-^)') 

as X ^ 00. 

Remark 2.9. By taking the Laplace transform of (2.17), one can obtain an asymptotic 
expansion of the (regularised) heat trace as t ^ 0. However, it seems that using the 
approach of [1] and [2], it is possible to obtain even stronger results (the pointwise 
asymptotic expansion of the heat kernel) . 
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Remark 2.10. Of course, formula (3.1) cannot be differentiated; moreover, we do not 
even know if in the almost periodic case N{X) is strictly increasing. However, in the 
periodic Schrodingcr case there are some results on the high-energy behaviour of the 
(non-integrated) density of states, see e. g. [11]. 

Given Conditions B and C, we want to introduce the following definition. We say 
that a non-negative function / = /(p) = /(p; k, k) satisfies the estimate /(p) ^ p^"*" 
(resp. /(p) ^ p*^^), if for each positive e and for each k we can achieve /(p) ^ p^ (resp. 
/(p) ^ p~^) for sufficiently large p by choosing parameter k from Conditions B and C 
sufficiently large. For example, we have 

(2.19) R{p) ^ pO+, 

card© ^ p^"*", s(p) ^ p°~, and r(p) ^ p^^ . 

Throughout the paper, we always assume that the value of k is chosen sufficiently 
large so that all inequalities of the form p^^ ^ p^ or p°~ ^ p^^ we encounter in the proof 
are satisfied. 

The next statement proved in [5] is an example of how this new notation is used. 

Lemma 2.11. Suppose, 0, n^, . . . , fj,^ G 0^, the set {^j} ^■^ linearly independent, and 
= Yl'j=ibjfJ'j- Then each non-zero coefficient bj satisfies 

p'-^\bj\^p'^. 

In this paper, by C or c we denote positive constants, the exact value of which can 

be different each time they occur in the text, possibly even in the same formula. On 
the other hand, the constants which are labeled (like Ci, C3, etc) have their values being 
fixed throughout the text. Given two positive functions / and g, we say that f ^ g, oi 
9^ f:Or g = 0{f) if the ratio g/f is bounded. We say / x 51 if / > 51 and / < g. 

We will also need a number of auxiliary constants. Let us choose numbers {aj}'j^i, (3, 

and <; satisfying 

(2.20) max{l - w + x/2, 1/2} <^<Q;i<Q;2<---<Q;d<i?<^<l 
(recaU (2.5)), and set 

(2.21) a := x/^. 

3. Reduction to a finite interval of spectral parameter 

To begin with, we choose sufficiently large po > Co (to be fixed later on) and for 
n G N put pn '■= 2p„_i = 2"po. We also define the intervals := [pn,4:pn]. The proof of 
Theorem 2.6 will be based on the following lemma: 
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Lemma 3.1. For each M G M there exist L > and a finite subset Jq G J such that for 
every n & N and p & 

d 1 L [^±M] 

(3.1) 7V(/-) = J2J2 H E ^'.^'in, M)p'^+(^-^-)'^+-+-+'-^- In^ p + 0(p;^). 

9=0 h=0 ti,...,theJo j=o 

Here, C^g^f^j'' {n, M) are some real numbers satisfying 

(3.2) C;-;''(n,M) = 0(p;^/^''+«). 

The constants in the O-terms do not depend on n (but they may depend on M). 

Remark 3.2. Note that (3.1) is not a 'proper' asymptotic formula, since the coefficients 
are allowed to grow with n (and, therefore, with p). 

Some of the powers of p on the right hand side of (3.1) may coincide. In order to avoid 
the ambiguity let us redefine coefficients C^'^'^'''' (n, M) in such a way that, for any given 
values of q and d + {2 — 2w)h + ii + ■ ■ ■ + ih — j , only the coefficient with the minimal 
possible value of h and maximal possible values of j, ci, . . . ,ih (in this order) is nonzero. 
Note that these new coefficients still satisfy (3.2). 

Let us prove Theorem 2.6 assuming that we have proved Lemma 3.1. Let M be fixed. 
Denote the sum on the right hand side of (3.1) by Nn{p'^^)- Then, for n ^ 1, whenever 
p e In-i n /„ = [pn, 2pn], wc havc: 

Nnip'n - Nn-i{p'n 

^^'^^ = E E E E ^^^r(^^ M)/+(^-^-)^+-+-+--^- In^ p + 0(p;^), 

q=0 h=0 ii,...,ih&Jo j=0 

where 

^ntin, M) qV/n^, M) - C-'/^n - 1, M). 
On the other hand, since for p e n /„ we have both A^(p2«') = A^„(p2"') + 0{p-^) 
and N{p^'") = A^„_i(p2«') + 0{p-^), this implies that 

d 1 L [ rf+M j 

(3-4) EE E E ^;vr(^,M)p<^+(2-2-)'^+-+-+^''-nnv=o(p;^). 

Claim 3.3. For each combination of indices present on the right hand side of (3.3) we 
have: 

(3.5) t'g'htin, M) = 0(p^-M-rf+(2«'-2)/^-M--.. 

Proof. Put y :— Pnj P and let 

(3.6) r-;'^(n,M) ■.= p'^^'^^^^-^^^^-^^^^-^^^-^Y.^^^ 
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Then by (3.4) for y e [1/2,1] 



d-1 L IT^I 



(3.7) P{y) -.= ^211 E E <.;''(^,M)|/^-'^+(2-2)'^----Mn^y = 0(1). 

p=0 h=0 (,i,...,(,hGJo j=0 

Let US denote by hi, . . . ,hT the functions yi-''+(2«'-2)/i-ti In''?/ entering the sum in 

(3.7) with non-zero coefficients; these functions are hnearly independent on the interval 
[1/2,1]. Therefore, there exist points £ [1/2,1] such that the determinant of 
the matrix {hj{yi))^- is non-zero. Now (3.7) and the Cramer's Rule imply that the 
values of Tp^h"j'''in', M) are fractions with a bounded expression in the numerator and a 
fixed non-zero number in the denominator. Therefore, 

(3.8) r-;''(n,M) = 0(l). 
Thus, choosing p = d — 1 in (3.6), we obtain 

td-ihji^^M) = o(p^-^-'^+(2«'-2)'^-'i--^'^). 

Now we can put p — d — 2 into (3.8) and obtain 

tt'ikM^M) = 0(p^-^-'^+M''----Mnp„). 

Continuing this process until p = 0, we obtain (3.5). □ 

Thus, for j < M + d + {2 -2w)h + Li + ■■■ + Lh, the series Y.Z=o ^qhT^'^^ ^) 
absolutely convergent; moreover, for such j we have: 

n 

C^j" in, M) = C;-/'' (0, M) + Y. ^^ht ^) 

ra=l 

oo 

=: C;V/''(M) + o(p^-^-'^+(2«'-2)'*-'i--^. In-'-^-Vn), 

where we have denoted C'^X't{M) := C;V/''(0, M) + Em=i CTl^^ Foi" bigger 
values of j we use (3.2) and (2.20) to obtain 

j- d+M j 

EE E E |C-;''(n,M)|p'^+(^-^-)'^+-+---+'-^-lnV 

9=0 /i=0 /.i,.--,the^o j>M+d+(2-2«;)/i+tiH ht/j 

^ E E E ln« p„ < In'^-^ p.. 

5=0 h=0 ti,...,(.^eJo 
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Thus, when p E In, we have: 

d-l L [M+d+{2-2w)h+ii+-+ih] 
qr=o h=0 ti,...,t;,eJo J=0 

Since the constants in O terms do not depend on n, it is sufficient to choose 

M:= [(<j(i + K)/(l-^)] +1 

to get (2.17) for all p ^ po- 

The rest of the paper is devoted to proving Lemma 3.1. The first step of the proof is 
fixing n and fixing large k and k. The precise value of k will be chosen later; the only 
restriction on it will be to satisfy inequality (9.9) (it says that the more asymptotic terms 

wc want to have in (3.1), the bigger k wc need to choose; note that the choice of k does 
not depend on k). We will have several requirements on how large k should be (most of 
them will be of the form p^ < or p°~ > p~^); each time we have such an inequality, 
we assume that k is chosen sufficiently large to satisfy it. 

Remark 3.4. Our choice of k will only depend on M, w, x, and the constants introduced 
in (2.20). The set Jo in Lemma 3.1 can be chosen to be 

(3.9) Jo := Jn [x-d-M- l,x]. 
The first requirement on k we have is that 

(3.10) k >d + M + x{d + M)/{w - x) -2w. 
After fixing k and k we get Rq from Condition B. Then, taking 

(3.11) Po > Ro 

and fixing n, we choose and J so that Conditions B and C are satisfied for p := 4p„. 
Without loss of generality we may assume that J D Jq. Then we introduce an auxiliary 
pseudo-differential operator B with the symbol b given by (2.8). 

Prom now on we prove Lemma 3.1 for S = S and with Jq replaced by J. However, 
in view of (2.13) and (2.20), the results with J and Jo arc equivalent. Afterwards, 
in Section 11 we will prove that the asymptotics (3.1) for the original B follows from 
Condition B and (3.1) for B. 

4. Pseudo-differential operators 

Most of the material in this and several subsequent sections is very similar to the cor- 
responding sections of [5] and [6], as are the proofs of most of the statements. Therefore, 
we will often omit the proofs, instead referring the reader to [5], [9], and [6]. 
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4.1. Classes of PDO's. Before we define the pseudo-differential operators (PDO's), 
we introduce the relevant classes of symbols. Let b — 6(x, x, ^ e R*^, be an almost- 
periodic (in x) complex-valued function and, moreover, for some countable set of 
frequencies (we always assume that is symmetric and contains 0; starting from the 
middle of this section, will be assumed to be finite) 

(4.1) 6(x,O = 5]K0,Oee(x), 

060 

where 

6(0,0 :=M.(6(x,0e_,(x)) 

are the Fourier coefficients of 6(-,^) (recall that M is the mean of an almost-periodic 
function). We always assume that (4.1) converges absolutely. Let us now define the 
classes of symbols we will consider and operators associated with them. For ^ e let 
(C) V^^+W- We notice that 

(4.2) (l + r7)<2(0(r7), V^,r/eR''. 

We say that a symbol b belongs to the class = Sa(/3) = Sa{(3, 0), if for any / ^ and 
any non-negative s e Z the conditions 

(4.3) l^lj"^^ :=maxy(6»)^sup(0(-"+l^l)^|D|6(6>,0| < oo, \s\ = s, + S2 + ■ ■ ■ + s^, 

6»g© 

are fulfilled. The quantities (4.3) define norms on the class Sq. Note that Sq is an 
increasing function of a, i.e. C for a < 7. 

Given e R'^, let us introduce a hnear map Ve on symbols which acts according to 
the rule 

(4.4) (y^)(cl>, a(cl>, ^ + 0)- a(cf>, 0- 
If the Fourier transform of the symbol is factorized, i.e. 

Q 

q=l 

then the action of Ve can be written as a sum of actions on each factor separately: 

Q q-1 Q 

(4.5) (v;^(0,o=En«^(*^'^+^)(^)('^'^) n «^('^'^)- 

q=l 1=1 s=q—l 

For later reference we mention here the following convenient bound that follows from 
definition (4.3) and property (4.2): 

(4.6) J2{^ysup{0^--+'^'^^[\BI(VJ){9,^)\)^C\^^^^^ s=\sl 
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with a constant C depending only on a, s, and /3. The estimate (4.6) implies that for all 
T7 with 1 77 1 ^ C we have a uniform bound 

Now we define the PDO Op(6) in the usual way: 

(4.7) Op(6)«(x) = (27r)-'^/2y" Kx,0e^^"(3^H)(0c?€, 

the integral being over W^. Under the condition 6 G the integral on the r.h.s. is clearly 
finite for any u from the Schwarz class S(]R''). Moreover, the property 6 G So guarantees 
the boundedness of Op(6) in L2(R''), see Proposition 4.1. Unless otherwise stated, from 
now on S(M'^) is taken as a natural domain for all PDO's when they act in L2(M'^). 

Applying the standard regularization procedures to definition (4.7) (see, e.g., [7]), we 
can also consider the action of Op (6) on the exponentials e^,, u ^ W^. Namely, we have 

Op(6)e^ = ^ h{d, i')e^+g. 

This action can be extended by linearity to all quasi-periodic functions (i.e. finite linear 
combinations of with different u). By taking the closure, we can extend this action 
of Op(6) to the Besicovitch space B2(M'^). This is the space of all formal sums 

00 00 
^^aje0^.(x), with |ajp < +00. 

3=1 j=i 

It is known (sec [7]) that the spectra of Op(6) acting in L2(M'^) and B2(M'^) are the same, 
although the types of the spectra can be entirely different. It is very convenient, when 
working with the gauge transform constructions, to assume that all the operators involved 
act in B2(M'^), although in the end we will return to operators acting in L2(]R'^). This 
trick (working with operators acting in 82 (K"')) is similar to working with fibre operators 
in the periodic case in the sense that we can freely consider the action of an operator on 
one, or finitely many, exponentials (2.3), despite the fact that these exponentials do not 
belong to our original function space. 

Moreover, if the order a = then by continuity this action can be extended to all of 
B2(M'^), and the extension has the same norm as Op(6) acting in L2 (see [7]). Thus, in 
what follows, when we speak about a pseudo-differential operator with almost-periodic 
symbol acting in 82, we mean that its domain is either whole 82 (when the order is 
non-positive), or the space of all quasi-periodic functions (for operators with positive 
order). And, when we make a statement about the norm of a pseudo-differential operator 
with almost-periodic symbol, we will not specify whether the operator acts in 12(1^^^) or 
82 (R'^), since these norms are the same. 
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4.2. Some basic results on the calculus of almost-periodic PDO's. We begin by 
listing some elementary results for almost-periodic PDO's. The proofs are very similar 
(with obvious changes) to the proof of analogous statements in [9] . 

Proposition 4.1. Suppose that h satisfies (4.1) and that |&|o°o < oo. Then Op (6) is 
hounded in both Ui^'^) and B2{W^) and \\ Op(6)|| ^ 

In what follows, if we need to calculate a product of two (or more) operators with some 
symbols bj G Saj{@j) we will always consider that bj G Sa^C^j &j) where, of course, all 
extra terms are assumed to have zero coefficients in front of them. 

Since 0^{b)u G S{W^) for any 6 G and m G S(R'^), the product Op(&)Op(£/), 
b G Sa(0i),5f G S^(02), is well defined on S(]R'^). A straightforward calculation leads to 
the following formula for the symbol bo g oi the product Op(6) Op(5f): 

(6o^)(x,0 = E + 

060i,<^66)2 

and hence 



(4.8) {b7g){x, = E ^(^' ^ + '^)^(^' 0, X e 01 + 02, e e 

e+<t>=x 



We have 

Proposition 4.2. Let b G Sa(©i), g G S^(02). Then bo g E Sa+^(©i + ©2) and 

with the constant C depending only on I, a, and s. 

We are also interested in the estimates for symbols of commutators. For PDO's 
A,-^i, I ^1,2,...,N, denote 

ad(A; ^'i, ^2, • • • , ^jv) := i [ad{A; ^1, ^2, • • • , ^iv-i), ^jv] , 
ad{A;'^) ■=i[A,% ad^(A; ^) := ad(A; ^, ^, . . . , ^), ad°(A; ^) := A. 

For the sake of convenience, we use the notation ad(a; ipi, ip2, ■ ■ ■ , V'at) and ad^(a, ip) for 
the symbols of multiple commutators. 
Let 

suppfe := {6> G R'' : b{9,-) ^ O}. 

It follows from (4.8) that the Fourier coefficients of the symbol ad{b,g) are given by 
(4.9) 

06(supp S)U(x-suppff) 
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Proposition 4.3. Let b G S„(0) and gj G 8^^.(0^), j = 1,2,...,N. Then 
ad(6; gi, . . . , gN) e S^(0 + 0^) with 

N 

7 = q; + ^(7^. -1), 

and 

N 

I (7.) 

\p,s+N-j+lJ 

where C and p depend on I, s, N, a and jj. 



\^^d{b;g^,...,g^)\l^^^C\b\^^l^ll\g,\ 



5. Resonant regions 

We now define resonant regions and mention some of their properties. This material 
is essentially identical to Section 5 of [5] , where the reader can find the proofs of all the 
statements of this section. ^ 

Recall the definition of the set = as well as of the quasi-lattice subspaces from 
Section 2. As before, by 0^ we denote the algebraic sum of k copies of 0; remember 
that we consider k fixed. We also put 0'^ := 0^ \ {0}. For each QJ e V we put 
— 1^ e QJ, = l|. For each non-zero e R*^ we put n(0) 

Let QJ G Vm- We say that is a flag generated by 03, if 5^ is a sequence QJj G Vj 
(j = 0, 1, . . . , m) such that 23j-i C 53^ and 2J„j = 03. We say that {i^jj^i is a sequence 
generated by ^ if Uj G Q3j093j_i and = 1 (obviously, this condition determines each 
up to multiphcation by —1). We denote by 3^(03) the collection of all fiags generated 
by 03. We put 

(5.1) Lj-.^p^^, 
recall (2.20). 

Let G 0'^. The resonant region generated by 6 is defined as 

(5.2) A(0):={|GM^ |(tn(0))| <Li}. 

Suppose, ^ G 3'(Q3) is a fiag and is a sequence generated by ^. We define 

(5.3) A(^):={£eR^ \{Cuj)\^Lj}. 

If dim 03 = 1, definition (5.3) is reduced to (5.2). Obviously, if C ^2, then A{^2) C 
A(^i). 

Suppose, 03 G Vj . We denote 

Si(Q3) := Uje:,(2j)A(5^). 
Note that Hi(3e) = R'^ and Hi (03) = A(0) if 03 G Vi is spanned by 6. Finally, we put 

(5.4) H(03) := Hi(03) \ ( Ukd^j 3i(il)) = 3^(03) \ ( Uud^i U^emH^))- 
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We call H(2J) the resonance region generated by QJ. Very often, the region H(X) is called 
the non-resonance region. We, however, will omit using this terminology since we will 
treat all regions S(5J) in the same way. 

The first set of properties follows immediately from the definitions. 

Lemma 5.1. (i) We have 

U5,evH(2J) = R''. 

(a) $, G Hi(QJ) ijf G f^(2J), where f2(5J) C 2J zs a certain bounded set (more 
precisely, = Hi(Q3) n 5J C 'B{mL^) z/dim2J = m). 

(Hi) Hi(M^) = H(M'^) IS a hounded set, H(M'^) C 'B>{dLd); all other sets Hi(QJ) are 
unbounded. 

Now we move to slightly less obvious properties. From now on we always assume that 
Po (and thus p„) is sufficiently large. We also assume, as we always do, that the value of 
k is sufficiently large so that, for example, Ljp^ < I/j+i. 

Lemma 5.2. Let G V. Then (3i(QJ) n Hi(il)) c Hi(Qn), where OH := QJ + 11 
(algebraic sum). 

CoroUciry 5.3. (i) We can re-write definition (5.4) like this: 

3(2J) :=3i(QJ)\(Uu522, Hi(il)). 

(ii) If^y^iX, then S(23) n H(il) = 0. 

(Hi) We have = U(j7evS(2J) (the disjoint union). 

Lemma 5.4. Let G Vm and C 2U G Vm+i- Let /x be (any) unit vector from 
2IT0QJ. Then, for ^ G Si(QJ), we have ^ G Si(2n) if and only if the estimate ~ 
^ Lm+i holds. 

Lemma 5.5. We have 

Hi(QJ) n UiiD<»Hi(il) = Hi(QJ) n U2HD2,, dim2D=l+dim2l3i(2IJ). 

Corollary 5.6. We can re-write (5.4) as 

(5.5) H(QJ) := Hi(QJ) \ ( U2BD2J,dim2U=l+dim2J Si(2IJ)). 

Lemma 5.7. Le^ QJ G V anc? G 0j^. Suppose that ^ G S(QJ) and both points ^ and 
i + e are inside A(0). Then 6 e ^ and ^ + 6 e H(QJ). 

Definition 5.8. Let 6,61,62,. ■■ ,Oihe some vectors from 0'^, which are not necessarily 
distinct. 

(a) We say that two vectors ^, 77 G M'^ are 6-resonant congruent if both ^ and rj are 
inside A{6) and — t]) = 16 with / G Z. In this case we write ^ O r; mod 6. 

(b) For each ^ G M*^ we denote by T0(^) the set of all points which are 0- resonant 
congruent to ^. For 7^ we say that Tei^) = if £ ^ A{6). 
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(c) We say that $, and r] are 61,62, ■■■ , 6i-resonant congruent, if there exists a se- 
quence e'R.'^,j^O,l,...,l such that £0 = ^, = rj, and E Te^{^j_i) for 
i = l,2,...J. 

(d) We say that 77 G M'' and ^ G M'^ are resonant congruent, if either ^ = t] or $, and 
77 are 61,62, ... , ^/-I'esonant congruent with some 61,62, ... ,61 G ©'^. The set 
of all points, resonant congruent to ^, is denoted by T(^). For points rj E T(^) 
(note that this condition is equivalent to ^ G T{rj)) we write T7 ■<->■ ^. 

Note that T(^) = {^} for any ^ G H(X). Now Lemma 5.7 immediately imphes 

Corollary 5.9. For each ^ G H(QJ) we have T(^) C H(QJ) and thus 

H(Q3) = U^eH(^u)T(0. 

Lemma 5.10. The diameter o/T(£) is bounded above by mLm, if ^ E S(QJ), QJ G V^. 

Lemma 5.11. For each ^ G H(2J), 5J 7^ R*^, ^/le set T(£) is finite, and cardT(^) is 
bounded uniformly m £ G M*^ \ H(]R'^). 

6. Description of the approach 

We first prove (3.1) assuming that the symbol 6 of i? is replaced by b which satisfies 
(2.8). In particular, it belongs to the class Sq,. At the end, in Section 11, we will use 
(2.14) to show that Theorem 2.6 holds as stated. 

For any set C C by CP(C) wc denote the orthogonal projection onto span{e|}|ge in 
B2(M'^) and by IP^(C) the same projection considered in L2(M'^), i.e. 

(6.1) 7^{e) ^TleS", 

where 5" is the Fourier transform and le is the operator of multiplication by the indicator 
function of C. Obviously, T^(C) is a well-defined (respectively, non-zero) projection iff 
6 is measurable (respectively, has non-zero measure). Let us fix sufficiently large n, and 
denote (recall that A„ = pf^) 

(6.2) X„ := G R', lil'"" G [(5/6)2- A„, A„]}. 
We also put 

Lemma 5.10 implies that, if po is big enough, 

(6.3) for each ^ G ^1 we have l^l^"' G [(2/3)2^A„, 62^A„] . 
In particular, we have 

(6.4) A n H(R'^) = 0. 
Let us define 

A:^{^^A, 1^1'" <A4 
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and 

(6.5) A:={^^A, ||P"'>A„}. 

We now plan to apply the gauge transform as in Sections 8 and 9 of [5] to the operator 
H. The details of this procedure will be explained in Sections 8 and 9; here, we just 
mention that we are going to introduce two operators: Hi and H2. The operator Hi 
is unitary equivalent to H: Hi = U^^HU, where U = e** with a bounded pseudo- 
differential operator with almost-periodic coefficients (then Lemma 2.3 implies that 
the densities of states of H and Hi are the same). Moreover, Hi = H2 + R^-, where 

(6.6) l|i?dl<Pn''+" 

and H2 = (—A)"" -|- is a self-adjoint pseudo-differential operator with symbol l^p"' + 
w^(x, ^) which satisfies the following property: 

(6.7) w^(6>,0 = 0, if (|^A(6>) k^eA), or + 6> ^ A(6>) & | e ^), or (6> ^ 0^). 

Wc can now use a simple statement which follows from Lemma 2.3 and Remark 2.7: 

Lemma 6.1. Suppose, Hi and H2 are two elliptic self-adjoint pseudo-differential oper- 
ators with almost-periodic coefficients such that \\Hi — H2\\ < p~'^^'^'^~'^ . Suppose that 
N{H2;p^'^) satisfies asymptotic expansion (3.1). Then N{Hi;p^'^) also satisfies (3.1) 
with the same coefficients. 

This means that it is enough to establish the asymptotic expansion (3.1) for the 

operator H2 instead of H. Condition (6.7) implies that for each ^ E A the subspace 
J'(T(^)) B2(M'^) is an invariant subspace of H2; its dimension is finite by Lemma 5.11. 
We put 

-f^2(^) -f^2|?(T($))B2(m'*)- 

Note that the subspaces 7{A)B2{R'^) and ?{A)B2{R'^) are invariant as well; by H2{A) and 
H2{A) we denote the restrictions of H2 to these subspaces; we also denote by H2{A) the 
restriction of H2 to 7{A)B2{W^). If wc consider the operator H2 acting in L2(]R'^), then 
?^{A)l2{M.'^), 7^{A)L20^'^), and {A)l2{M.'^) are still invariant subspaces. It follows 
from (6.2) - (6.5) that UH2{A)U* < {^/Qf'"\J and UH2{A)U* > ^'^'"XJ. 

For each E A the operator i^2(^) is a finite- dimensional self-adjoint operator, so 
its spectrum is purely discrete; we denote its eigenvalues (counting multiplicities) by 
-^1(0 ^ '^2(^) ^ . . . ^ AcardT(^)(^). Ncxt, WC list all poiuts 77 G T(^) in increasing order 
of their absolute values; thus, we have put into correspondence to each point 77 G T(^) a 
natural number t = t{r]) so that t{rj) < t{ri') if |r7| < |r7'|. If two points r] = {rji, . . . ,rjd) 
and T]' = {r]'i, . . . ,r)'^) have the same absolute values, we put them in the lexicographic 
order of their coordinates, i.e. we say that t{r}) < t{ri') if r]i < rj'i, or r)i = ri[ and 
V2 < ''721 stc. Now we define the map g : A ^ which to each point r] E A brings 
into correspondence the number At(^) (T(r7)) . This map is an injection from A onto the 
set of eigenvalues of H2, counting multiplicities (recall that we consider the operator 
H2 acting in B2(R''), so there is nothing miraculous about its spectrum consisting of 
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eigenvalues and their limit points). Moreover, all eigenvalues of H2 inside the interval 
[(7/8)2^ A„, (9/2)2^ A„] have a pre-image under g. We define 

(6.8) ^(^) := l^p^-, for ^ e R'* \ A 

Arguments similar to the ones used in [6] show that is a measurable function. 
We introduce 

Lemma 6.2. For X e [A„,4^"'A„] being a continuity point of N{\;H2) we have: 

(6.9) A^(A;i/2) = (27r)-'^voigA. 

Since points of continuity of N{X) are dense, the asymptotic expansion proven for such 
A can be extended to all A G [A„,4^'"A„] by taking the limit. Thus, our next task is to 
compute volS;^. Let us put 

A+{p) := e R^ g{^) < p'^ < 1^^} 

and 

A-(p) := {^eR^ </"'<^(0}. 

Lemma 6.3. 

(6.10) voKSa) = + vol^+(p) - yol A-{p), 
where Ud is the volume of the unit ball in R^. 

Proof We obviously have Sa = 'B{p)UA+{p) \A-{p). Since A'ip) C S(p) and A+{p) n 
!B(p) = 0, this implies (6.10). □ 

Remark 6.4. Properties of the mapping g imply that A^{p) UA^lp) C A. Thus, in 
order to compute A^(A), we need to analyze the behavior of g only inside A. 

We will compute volumes of A^{p) by means of integrating their characteristic func- 
tions in a specially chosen set of coordinates. The next section is devoted to introducing 
these coordinates. 

7. Coordinates 

In this section, we do some preparatory work before computing voiyi''=(p). Namely, 
we are going to introduce a convenient set of coordinates in S(5J). Let 2J G be fixed; 
since A^{p) n3(R'^) = 0, we will assume that m < d. Then, as we have seen, ^ G Hi(QJ) 
if and only if G Q(QJ). Let {11^} be a collection of all subspaces iXj G V^+i such that 
each iij contains QJ. Let = /ij(5J) be (any) unit vector from iij QJ. Then it follows 
from Lemma 5.4 that for ^ G 3i(QJ), we have ^ G Hi(ilj) if and only if the estimate 
|(€)A*j)| = |(€aj^)Mj)| ^ Lm+i holds. Thus, formula (5.5) implies that 

S(QJ) = G R^ G Q(QJ) & Vj |(^^x,/u,.(QJ))| > L^+i}. 
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The collection |yLij(Q3) j obviously coincides with 

{n(0^x), ee&-,\^}. 

The set 3(QJ) is, in general, disconnected; it consists of several connected components 
which we will denote by |S(QJ)p|^__^. Let us fix a connected component S(53)p. Then 

for some vectors C {±Atj} we have 

we assume that {/i^ (j9)}^^^ is the minimal set with this property, so that each hyperplane 

has a non-empty intersection with the boundary of H(QJ)p. It is not hard to see that 
Jp ^ d — m. Indeed, otherwise 3(QJ)p would have non-empty intersection with Hi(QJ') 
for some QJ', QJ C QJ'. We also introduce 

Note that our assumption that S(QJ)p is a connected component of 3(2J) implies that 
for any ^ G S(2J)p and any G 0^ \ QJ we have 

We also put 

K := d — m — 1. 

Without loss of generality we may (and will) assume that the number Jp of 'defining 
planes' is the minimal possible, i.e. Jp — K + 1. Indeed, the argument presented in 
Section 11 of [5] explains how to derive the result for arbitrary S(QJ)p, assuming we have 
proved it in the case Jp = K + 1. 

If Jp = K + 1, then the set is linearly independent. Let a = a(p) be 

a unique point from QJ-*- satisfying the following conditions: {a, p,j{p)) = Lm+i, j = 
1, . . . , K + 1. Then, since the determinant of the Gram matrix of vectors is > 

by (2.16), we have 

We introduce shifted cylindrical coordinates in H(QJ)p. These coordinates will be denoted 
by ^ = (r; X). Here, X = (Xi, . . . , Xm) is an arbitrary set of cartesian coordinates 
in f2(QJ). These coordinates do not depend on the choice of the connected component 
H(QJ)p. The rest of the coordinates (r, #) are shifted spherical coordinates in QJ-*-, cen- 
tered at a. This means that 

r(^) = - a| 

and 

$ = n(^2jx - a) e ^21^. 
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More precisely, $ G Mp, where Mp := {n(^2ix ~ ^i), ^ G H(Q3)p} C ^^jx is a i^T- 
dimensional spherical simplex with K + 1 sides. Note that 

Mp = {n(^5jx - a), ^ G H(Q3)p} = {n(^2jx - a), Vj (^^J^, > Wi} 

= {n(r7), r7 := l^jx - a G QJ^, Vj (77, > O} = 5^x n H(9J)p. 

We will denote by d$ the spherical Lebesgue measure on Mp. For each non-zero vector 
fjt G QJ"*", we denote 

W(/u) := {r7 G QJ^ (r7,iu) = 0}. 

Thus, the sides of the simplex Mp are intersections of W[fj,j{p)) with the sphere -S'^jx. 
Each vertex v = Vt,t = 1,.. .,7^ + 1 of Mp is an intersection of S'^jx with K hyperplanes 
W(/ij(p)), j = 1, . . . , K + 1, j ^ t. This means that is a unit vector from QJ-*- which 
is orthogonal to j = 1, . . . , K + 1, j ^ t; this defines v up to a multiplication 

by-1. 

Lemma 7.1. Let iii and ii2 be two strongly distinct subspaces each of which is a linear 
combination of some of the vectors from {//^ (p)}. Then the angle between them is not 
smaller than s{pn)- In particular, all non-zero angles between two sides of any dimensions 
of Mp as well as all the distances between two vertexes and Vt, t ^ t, are bounded 
below by s{pn)- 

Lemma 7.2. Let p be fixed. Suppose, G 0^ \ QJ and 0^± — Ylf=i bjfljip). Then 
either all coefficients bj are non-positive, or all of them are non-negative. 

By taking sufficiently large k wc can assure that the diameter of Mp does not exceed 
(lOOci^)^^ We put $g := f - 0(|2jx - a, q = 1,. ..,K + 1. The geometrical 

meaning of these coordinates is simple: $q is the spherical distance between $ = n(^(jjx — 
a) and W(/ig(p)). The reason why we have introduced $5 is that in these coordinates 
some important objects will be especially simple (see e.g. Lemma 7.5 below) which 
is very convenient for integration. At the same time, the set of coordinates (r, {^g}) 
contains K + 2 variables, whereas we only need K + 1 coordinates in QJ^. Thus, we 
have one constraint for variables $j. Namely, let {ej}, j = 1, . . . + 1 be a fixed 
orthonormal basis in QJ-*- chosen in such a way that the K + 1-st axis is directed along 
a, and thus passes through Mp. Then we have e^ = X^^t^ with some matrix {%;}, 
j, Z = 1, . . . , X + 1, and jli — iii{p). Therefore (recall that we denote rf :— — a), 

K+l 

Vj = iv, ej) = r ^ ajg sin 

9=1 

and, since r^(^) = jryp = X^jLV^I' ^^^^ implies that 

K+l K+i 2 

J2 {J^^^i^^'^^i) = 1' 
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which is our constraint. 
Let us also put 

(7.2) = 

Then we can write the surface element in the coordinates {77^} as 

dri[ . . . dri'j^ drj'^ . . . drj'j^ 



d^ 



^^+1 (i-Ef=i(^;)^) 



where the denominator is bounded below by 1/2 by our choice of the basis {ej}. It 
follows from our choice of the coordinates and (7.2) that 

K+l 

(7.3) (a, $) = (a, 11(77)) = \^Wk+i = |a| Yl (^K+igSin 

q=l 

Lemma 7.3. For each p, I we have \api\ ^ s{pn)~^- 
Lemma 7.4. We have max^ sin $^ (77) ^ s{pn)d~^^'^. 

The next lemma describes the dependence on r of all possible inner products {^,0), 

Lemma 7.5. Let ^ e H(QJ)p, QJ G and e 
(i) If E^, then {$,, 0) does not depend on r. 
(a) If and 0<x3± = Tliq^qi^qip)} then 

it 0) = (X, 6>2j) + Lm+i Y.^q + r{i) bq sin 

g g 

In the case (ii) all the coefficients bg are either non-positive or non-negative and each 
non-zero coefficient bg satisfies 

(7.4) pI < \bg\ < p^- 

8. Partition of the perturbation 

The symbols we are going to construct in this section will depend on p„; this depen- 
dence will usually be omitted from the notation. 
Let w e C°°(M) be such that 

1, ^^1; 

0, 21/20. 



^.1) O^zu ^1, w{z) 
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:= zu 




:= 1 - 




:= 1 - 



and 
(8.3) 



Ce(0 



1(^,^ + 0/2)1 



Pn|6l| 



MO := 1-^(1). 



Remark 8.1. Note tliat eg + £q + ig = 1. Tlie function is supported on the set 
1^ + 0/2\ ^ 7pn, and £g is supported on tlie set |^ + d/2\ ^ Pn/2. Tlie function ee is 
supported in the shell pn/S ^ |^ + 0/2| ^ Spn- 

Using the notation £g for any of the functions ig or ig, we point out that 

MO = ^-0{^ + e), CeiO^C-oi^ + d). 
Note that the above functions satisfy the estimates 



(8.4) 



|D|e,(0| + |D|£,(0| 



|D|^.(^)| + |D|C.(€)|<Pn''l^'. 

Now for any symbol b G Sa(/3) we introduce five new symbols: 

&^'(x,£;Pn) 



&^^(x,£;Pn) 
^'''(x,^;p„) 

&''(x,^;Pn) 
&''(x,^;Pn) 



?e0' 



6»e0' 

6''(ePn):= 6(0,0- 

The superscripts here are chosen to mean, respectively: 'large energy", 'non-resonant', 
'resonant', 'small energy' and 0-th Fourier coefficient. The corresponding operators are 
denoted by B^^, B^^, B^, B^^, and B°. By definitions (8.1), (8.2) and (8.3) 

(8.5) b^b" + b^^ + b^ + b^^ + b'^^. 
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The role of each of these operators is easy to explain. Note that on the support of the 
functions ¥^^{6, •; p„) and b^{0, •; p„) we have (using (2.19)) 

p./3-0(pr)^|^|^8p„ + 0(pr). 

On the support of b^^{6, •; pn) we have 

(8.6) |C|^p„/2 + 0(p°+). 
On the support of b^^{0, ■ ; p„) we have 

(8.7) m>7pn-0{p'i+). 

The introduced symbols play a central role in the proof of Lemma 3.1. As wc have 
seen in Section 6, due to (8.6) and (8.7) the symbols b^^ and b^^ make only a negligible 
contribution to the spectrum of the operator H near A = p^"' for p e In- The only 
significant components of b are the symbols b''^^, b^ and b°. The symbol b° will remain 
as it is, and the symbol b^'^ will be transformed in the next section to another symbol, 
independent of x. 

Under the condition b e Sa(/3) the above symbols belong to the same class Sq(/3) and 
the following bounds hold: 

i^^i!? + i^^'iff + + + \b''\ts ^ \b\ts- 

If b symmetric, then so are the symbols on the right hand side of (8.5). 

Let us mention some other elementary properties of the introduced operators. In the 
lemma below we use the projection J'(C), 6 C M*^ which was defined in Section 6. 

Lemma 8.2. Let b e Sa(/3) with some a e R. Then: 

(i) The operator B^^ is bounded and 

Moreover, 

(l - y(S(2p„/3)))s^^ = - y(S(2p„/3))) = 0. 

(ii) The operator satisfies the relations 

T(a3(p„/6))5* = 5*y(a3(p„/6)) = (/ - T(s(9p„)))£^ = - y(9p„))) = o, 

and similar relations hold for the operator B"^"^ as well. Moreover, ft'^^, b^ e S-y 
for any 7 e M, and for all I and s 

with the implied constant independent ofb andn ^ 1. In particular, the operators 
B-^^, B^ are bounded and 



(a) 
■ 



ASYMPTOTIC EXPANSION OF IDS 25 

(iii) 

9. Operators Hi and H2 

9.1. Preparation. As mentioned at the end of Section 3, we assume that the symbol b 
of B satisfies (2.8), and thus belongs to the class Sa{P) with a defined in (2.21). Our 
strategy is to find a unitary operator which reduces H = Hq + B, Hq :— (—A)"', to 
another PDO, whose symbol, essentially, depends only on ^. More precisely, we want to 
find operators and H2 with the properties discussed in Section 6. 

Repeating the calculations of Subsection 9.1 of [5] we find that H is unitarily equivalent 

to 

(9.1) H^ = Ho + yS"^ + If + F-''''^' + R-^ , 
where 

k k 

(9.2) n:=E^^+E^'' 

1=1 1=2 

B, := Op(6), 

i-i ^ 

(9.3) ^(-El E ad(0p(6);*,„*fe„...,*fe,), ^^2, 

i=l ki+k2+-+kj=l-l 

' 1 

(9.4) ^;-E-T E ad(i/o;*iki,*te,---,*fc,), i^2, 

j=2 ki+k2+-+kj=l 

Rj,:^ I dti I dt2--- /" %xp(-it^)ad^+^(i/;^)exp(it^)dt 
Jo Jo Jo 

k 



+ E ad(ii';*fe,,*fc„...,*fc.). 



3- 

i-l ki+k2+-+kj^k+l, 
kq^jk, q=l,-,j 

k 

p=i 

The symbols ijjj of PDO are found from the following system of commutator equa- 
tions: 

(9.5) ad(i7o;*i) + 5f^ = 0, 

(9.6) ad(//o; *;) + + 7]^^ = 0, 1^2. 
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By Lemma 8.2(ii), the operators Bf^^, Tf^''^ are bounded. This, in view of (9.5) and (9.6), 
imphes boudedness of the commutators ad(ifo;^i)> ^ ^ 1- Below we denote by yj^ the 
symbol of the PDO Y~^. 

9.2. Commutator equations. Put 

(9.7) U^) := eeiOMOil^ + - l^r)'' 

when 6 0, and Xo(€) '■— 0- 

We have 

Lemma 9.1. Let A = Op(a) be a symmetric PDO with a G S^,. Then the PDO ^ with 

the Fourier coefficients of the symbol ip{-x,^) given by 

(9.8) i^{e,^):^id{e,^)x0{^) 

solves the equation 

ad(//o; *) + Op(a^^) = 0. 

Moreover, the operator ^ is bounded and self-adjoint, its symbol ip belongs to S-y with 
any 7 e M and the following bound holds: 

The proof of this lemma is analogous to that of Lemma 4.1 of [6] and is based on the 
estimate 

1^ + 6>|2^ - = 1^1'"' ((1 + |^r'(2^ + 0) ■ ey - 1) X p2"'-2|6> . + 0/2)1 

which holds for ^ in the support of cetpe. 

Using Propositions 4.1, 4.2, 4.3, Lemma 9.1, and repeating arguments from the proof 
of Lemma 4.2 from [6] (with Oj :— j{a — 2 — {2w — 2)/3~^) + 1), we obtain the following 

Lemma 9.2. Let b e Sa{(3) be a symmetric symbol. Suppose that k is large enough so 
that r[pn)~^ < < pZ^^ ^ ^ and k satisfies 

(9.9) k > 2(M + ai3 + d- 2w)/{2w + 2/3 - a/3 - 2). 

Then ipj, bj, tj e S-y(/3) for any 7 e M and there exists sufficiently large po, such that 

(9.10) im^Pn''^''"-'- 

Remark 9.3. Note that the expression in the denominator of (9.9) is positive by (2.20) 
and (2.21). 

Now Lemmas 6.1 and 9.2 imply that the contribution of Rj;. to the integrated density 
of states can be neglected. More precisely, let be the operator with symbol 

(9.11) ^^(x,0:=y^(x,0-Z/f(x,0, i-e. w~,{0,^) ^ - eei^M^))- 
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We introduce H2 := + Wj^. Then, by (9.1) and (9.10), \\Hi - ifall < p'^^'^'^''^ 

and, moreover, the symbol wj. satisfies (6.7). This means that all the constructions of 
Section 6 are valid, and all we need to do is to compute vol Sa- 

Until this point, the material in our paper was quite similar to the corresponding parts 
of [5]. From now on, the differences will be substantial. 



9.3. Computing the symbol of the operator after gauge transform. The follow- 
ing lemma provides us with a more explicit form of the symbol y-j^. 

Lemma 9.4. We have y^{d,^) = for 6 ^ 0^. Otherwise, 

s=i ej,9s+ie@ p=i £)^',<^^'e©;,+i !^i,-,i^2s+p>0q=i 

(9.12) (t>j,<t>,+i,<t>'je&s+i E^i=s 

s 

X (V^)(0,+i, ^ + </>,+i) n (V^6)(0,-, £ + (l>j){V'^xe',M + <t>'i)- 
Here for 1/ e N 

(9.13) J2 V°:=C(^'^)({0,0}), 
T7i,...,T7^e© 

and, for e M.^, the action ofVg on symbols of PDO is defined in (4.4), whereas for 
any function f on 

{^efm ■.= f{i + e)-f{i). 

The coefficients C(^'f)({0, 0}) and C^^;^].,n, {{6, 0}) depend on s, p and all vectors 0, 
0j, ds+i, <pj, <Ps+iJ ^'jj ^qj 'Pq (dnd ou r}i,...,r}^ if these subscripts are present). 

Moreover, these coefficients can differ for each particular , z/ G Nq- At the same 
time, they are uniformly bounded by a constant which depends on k only. We apply the 

convention that ng=i ^^'^^0q'^e'^){'i + 0g) = 1- 

Proof. We will prove the lemma by induction. Namely, let ^ ^ 2. We claim that: 
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1) For any m = 1, — 1, tpm{G, ^) = for ^ Otherwise, 

m p— 1 

E E E E n(^'^'^«;')(€+<) 

e^ee p=i 6»^',<^^'g0;„ j^i,...,j^2m+p-i>0g=i 
(9.14) e^ee; 



9'e©' 



2) For any s = 1, 1 and any ki, . . . ,kp {p ^ 1) such that ki -\- ■ ■ ■ -\- kp — 
ad ( Op(6)TSr, ■ ■ ■ , l^', = for e,+i. Otherwise, 

ad(Op(6)7ir,---,*O(^'0 

= E E E E n(v'^'^.^')(^+<) 

(9.15) </.^.,</>3+i,<;i^ee,+i is;gs;p-i E'^i=5 

X (V^6)(0,+i, £ + </>,+i) n (V^6)(0,-, £ + </',-)(V^^^0;.)(^ + </>;). 

3) For any s = 2, . . . ,i and any ki, . . . ,kp {p ^ 2) such that ki + • ■ ■ + kp — s, 
ad{Ho; ^Z^. . , *fej(6>, = for 6> ^ 0,. Otherwise, 

ad(i/o;*^..,*O(0,O 

= E E E E n(V'^'^.;')(c+<) 

ej,e,e& p=i e'^,<t>'^&&s vi,...,v2s+v-2^o q=i 

(9.16) <t>j,4>„<t>'je&s Kg^p-i E'^x=5-i 

X (V^)(0„ ^ + (/)J n (V^&)(0,-, £ + </>,)(V-^^X0j)(£ + 0;.). 

For £ = 2 assumptions l)-3) can be easily checked. Indeed, by (9.8), (4.9) and (4.5), 

ad ( O^; *i) (0, = E ikO.me-x.i-rx) ( V^^^) {$) 

xe©u(6>-©) 
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XG0U(0-0) 

- Mi)ee{mo, mo -x,i+ x)iyx^-xm 

- Mi)ee{mo.i)^){e - x.i)xB-x{i))- 

Now, we complete the induction in several steps. 

Step 1. First of all, notice that due to (9.3), (9.4), for any m = 2, . . . ,£ the symbol of 
Bm admits a representation of the form (9.15) with s = m — 1, and symbol of Tm admits 
a representation of the form (9.16) with s — m. Then it follows from Lemma 9.1 and 
(9.6) that '^i admits a representation of the form (9.14). 

Step 2. Proof of (9.15) with s^L Let /ci H ^kp^l. If p ^ 2. Then 

ad ( Op(6); , . . . , ^k,) = ad ad ( Op(6); ^k,,---, ; ^fc,) ■ 

Since ki + ■ ■ ■ + kp^i ^ £ — 1 and kp ^ £ — 1 wc can apply (9.14) and (9.15). Combined 
with (4.9) it gives a representation of the form (9.15). If p = 1 then ad ( Op(6); ^'f) 
satisfies (9.15) because of (4.9) and step 1. 

Step 3. Proof of (9.16) with s^£+l. Let A;i H ^kp^£ + l,p^2. U p ^ 3, then 

(cf. step 2) 

ad{Ho; ^'fc,, . . . , ^^fcj = ad ( ad(i/o; ^fei, • • ■ , *fcp_J; ^k,)- 

Since ki + ■ ■ ■ + kp-i ^ £, p — 1 ^2 and fcp ^ £ — 1 we can apply (9.14) and (9.16). 
Together with (4.9) it gives a representation of the form (9.16). If p = 2 then (see (9.6)) 

ad(i/o; = ad ( ad{Ho; ^k,) = - ad(i?2^ + T.f ; ^^J. 

Since ki ^ £ and /c2 ^ £, the representation of the form (9.16) follows from (4.9) and 
step 1. (Formally exceptional case ki — 1, k2 — £ can be treated separately in the same 
way using (9.5) instead of (9.6).) 
Induction is complete. 

Now, (9.15), (9.16) and (9.2), (9.3), (9.4) prove the lemma. □ 

10. Contribution from various resonant regions 

Let us fix a subspace QJ G V„j, ni < d, and a component Hp of the resonant region 3(QJ). 
Our aim is to compute the contribution to the density of states from each component 
Hp. Therefore, we define 

(10.1) ^+(p) :=^+(p)nHp and yi;(p) yi-(p) n H^ 

and try to compute 



(10.2) 



voiyi+(p)-voi^;(p). 
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Since formulas (6.10) and (6.4) obviously imply that 

d-l 

(10.3) voKSa) = ujap" + E E E(^ol-^p ('^) - vol^;(p)), 

m=0 QJeV^ p 

Lemma 3.1 would be proved if we manage to compute (10.2) (or at least prove that this 
expression admits a complete asymptotic expansion in p). 

Note that if ^ G Hp, then we also have that T(^) C Hp. We denote 

H^{i) := H2\h,, H^:=y(T(0)B2(M'') 

(recall that is an invariant subspace of H2 acting in B2(]R'^)). Suppose now that two 
points ^ and rj have the same coordinates X and $ and different coordinates r. Then 
^ e Hp implies e Hp and T(?7) = T(^) + (77— ^). This shows that two spaces and H^^ 
have the same dimension and, moreover, there is a natural isometry : — > given 
hy F : e,^ ^ e,^+(jj_^), u G T(^). This isometry allows us to 'compare' operators acting 
in and H^. Thus, abusing slightly our notation, we can assume that i^2(^) and -f^2(^) 
act in the same (finite dimensional) Hilbert space H(X, $). We will fix the values (X, ^>) 
and study how these operators depend on r. Thus, we denote by H2{r) = H2{r; X, $) 
the operator i?2(^) with ^ = (X,r, $), acting in H(X, $). 

Let W^{r) be the operator in H(X, $) with the symbol w^(x, £(X, r, ^)). According 
to formula (7.3), for any s ^ k — 1 and 6 e &s+i 

K+l 

(10.4) 1^ + = + 2r|a| E aK+iqSm + 2(^, 0) + |X|' + \a\^ + 

9=1 

This, together with (2.4), (2.10) and (2.11), implies that for |^ + </>| > Cq the coefficients 
^ + 0) can be represented as the absolutely convergent series 

00 K+l 

6(0, 1 + <^) ^ E E E ^i^X,. ^y~'^i *„)- , 

(10.5) ^gj 1=0 ni,...,nK+i>0 a=l 

niH hriK+l^' 

jiH Vjd^l 

where the coefficients satisfy 

In the next lemma, to facilitate the expansion of the RHS of (9.12) in a suitable form, 
we transform the denominator of Xe' (recall (9.7)). 

In the subsequent calculations we will use the generalized binomial coefficeints: 

1, J = 0; 

:i0.6) [^\:=l 1-^ 



.3) ' \-J[iP-^)^ jeN. 

fe=0 
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Lemma 10.1. For s ^ k — 1, <p' G 0^+1, 6' G Q^+i; and ^ in the support of e^'cpg' let 

^ j=2 ^-^^ ik=0 ^ ^ 9=1 

Then \D\ < p-i+"-+i+o+ and 



oo 

(10.7) (|£+(/)'+6>f"'-|£+<^f"')"' =«;-V-2«'(2(t6>')+2(<^',6»') + |6'f)"'^(-L')". 

a=0 

Proof. We introduce a shorthand 

K+l 

N := 2r|a| J] ai^+igSm$, + 2(e </>') + IX]" + |a|' + |0f . 

q=l 

Then by (generahzed) binomial formula and (10.4) we obtain 
1^ + </>' + 6>f^ - 1^ + 

= (l^r + 2(^, <f>' + e') + \ci>' + e'\Y - (i^r + 2(^, cf>') + icp'iY 

^{r' + N + 2{t e') + 2(0', 6>') + \e'\Y - (r' + N)^ 



(10.8) 



oo ^ X 

= r^- I ((iV + 2(£, d') + 2(</,', 0') + |6>f - 7V^) 



= wr^^-\2{t e') + 2(<^', e') + |6lf ) (1 + D). 

The estimate on \D\ follows from estimates (2.19) and (7.1), and Lemmas 7.3 and 5.1. 
Now (10.7) follows from (10.8). □ 

As we have seen from the previous sections, the symbol of the operator H2 satisfies 

(10.9) /i2(x, = III'" + «^^.(x, i) = {r' + 2r(a, + |a|2 + + ^/;^(x, 0, 

where wj^ are given by (9.11) and (9.12). 

Remark 10.2. In this section we assume that ^ G j^, so by (6.3) 2p„/3 ^ ||| ^ 6p„, 
and by Remark 8.1 all functions eg{^ + ■) from (9.11) and (9.12) are equal to 1. Note 
that if G 0^, G 0fc, and QJ, then (see Lemma 7.5 and (8.3)) ipei^ + </>) = 1. 
This means that all cut-off functions from (9.11) and (9.12) are equal to 1 unless G 5J. 
If, on the other hand, G 5J, then ipe{$, + 0) depends only on the projection and 
thus is a function only of the coordinates X. 

By Proposition 4.1, (9.11), Lemma 9.4, formulas (10.5) and (10.7), Lemma 7.5, and 
Remark 10.2, for r x p„ 
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This, together with (10.9), imphes 

Lemma 10.3. The operator H2{r) is monotonically increasing in r; in particular, all its 
eigenvalues \j(^H2{r)^ are increasing in r. 

Thus the function g(^^(X., r, $)) (defined in Section 6) is an increasing function of r if 
we fix the other coordinates of ^, so the equation 

has a unique solution for fixed values of X and we denote the r-coordinate of this 
solution by T — t{p) — t{p; X, so that 

(10.11) ^(^(X,T,$))=p^-. 

By To = 7o(p) — 'To{p', X, $) we denote the value of r for (—A)"', i.e. Tq is a unique 
solution of the equation 

||(X,To,*)|=p. 

Obviously, we can write down a precise analytic expression for Tq (and we have done this 
in [4] in the two-dimensional case) and show that it allows an expansion in powers of p 
and Inp, but we will not need it. The definition (10.1) of the sets Ap{p) implies that 
the intersection 

yi;(p)n{^(x,r,*), reR+} 

consists of points with r-coordinate belonging to the interval [to{p),t{p)] (where we 
assume the interval to be empty if tq > r). Similarly, the intersection 

/l;(p)n{^(X,r,$), rGM+} 
consists of points with r-coordinate belonging to the interval [r(p), To(p)] . Therefore, 

= = ^(X, r, #), X e Q(QJ), $ e Mp, r e [to(p; X, $), r(p; X, $)] } 

and 

•A;(p) = - €(X, r, *), X e 0(QJ), * e Mp, r e [t{p; X, *), ro(p; X, *)] }. 
This imphes that (recall that K — d — m — 1) 

volA+{p) -volA-{p) ^ / dX r^dr 

[ d^ [ dX{T{p; X, *)^+^ - ro(p; X, *)^+^) . 

Remark 10.4. Note that in the case K = the simplex Mp is degenerate and there is 
no integration in 
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Obviously, it is enough to compute the part of (10.12) containing r, since the second 
part (containing tq) can be computed analogously. We start by considering 

(10.13) / T(p;X,*)^+^dX. 

First of all, we notice that if ^, 77 e S(5J) are resonant congruent points then, according 
to Lemma 5.7, all vectors Oj from Definition 5.8 of equivalence belong to QJ. This 
naturally leads to the definition of equivalence for projections and ry^j. Namely, we 
say that two points u and /Lt from fl{^) are 2J-equivalent (and write u -H-^j /Lt) if v and 
H are equivalent in the sense of Definition 5.8 with an additional requirement that all 
6j G QJ. Then ^ •H- 77 implies ri^jj. For 1/ G il(2J) we denote by T(u(i/) the class 

of equivalence of u generated by -H-su. Then y<n{^<xs) ^ projection of T(^) to QJ and is, 
therefore, finite. 

Since T2j(i/) is a finite set for each u G i^(Q3), we can re- write (10.13) as 



(10.14) f r(p;X,*)^+i(iX= / (card Ya,(i/))"' V t(p;X,* 
and try to compute 



Remark 10.2, together with equations (10.9), (9.11), and (9.12), shows that H2{r) 
depends on r analytically, so we can and will consider the family H2{z) with complex 
values of the parameter z with 3fte 2; x p. Likewise, we analytically continue the function 
$(X,r,*) to 

(10.15) ^(X,z,*) := X + a + z$. 

We also introduce the analytic continuation | ■ |c of the modulus of vectors, so that 

(10.16) := + 2z{a, $) + jap + |Xp. 

Formulas (9.11) and (9.12) give matrix elements of i72(^) in an orthonormal basis even 
for complex z. 

We choose a contour 

(10.17) 7 := 1^ G C : =tp„ := (^8max{(2w-2)/3, 1}) Vnj 

to be a circle in the complex plane going in the positive direction. 

Estimates (10.10) remain valid after the analytic continuation: for all z inside and on 

7 

' ^ < o^-'+o+ 7 > 



(10.18) 



Lemma 10.5. For p E In = [Pn)4p„] all r(p;X, $) lie inside 7. These are the only 
zeros of the function det {H2{z) — p^^l) inside the contour. 
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Proof. Let r := 3fte z, y := '^m z. For y = the operator H2{r) is self-adjoint. Thus it 
has card 'T<u{iy) real eigenvalues. 

Now for r ^ p + tpn'^ {l + t/A)p relations (10.9), (7.1), Lemma 5.1(ii), and (10.18) 
imply 

H^ir) ^ (((l + V4)p)'"(l-0(p"-+^-^+o+)) -0(p-+°+))/. 

Thus by (2.20) and (2.21) for big p no eigenvalue of H2{r) can coincide with p^"'. 
Likewise for r ^ p — tpn ^ (1 — t/4)p for big p we have 



Ho(r) ^ 



(((l-t/4)p)^'"(l + 0(p"---^+°+))+0(p 



and no eigenvalue of H2{r) can coincide with p^"'. This implies that all the eigenvalues 
of H2{r) lie in the real interval (p — tpn,p + tpn)- By (10.11) and Lemma 10.3 these 
eigenvalues coincide with {r(p; X, $) : X e t^{u)Y 

It remains to show that H2{r + iy) is invertible for any nonzero y such that r + is 
inside or on 7. Relation (10.15), Lemma 5.1(ii), definition (5.1), and bound (7.1) imply 
that inside and on the contour 



and 



Hence 



£=(r + iy)(l + 0(p;i+"-+^+°+)) 
^-rg l^lc ^ (1 + 0(1)) arcsin(ip„/p) ^ (l + o(l)) arcsini ^ t{l + o(l)). 



\2w I 

Ic 



|2 I 

Ic 



p^'" and argl^l^' 



2|/(r + (a,*)) 



'^'^ — w arcsin , ' x yp ^, 



which implies that 
(10.19) 



> 



y\p 



2w-l 



Now for any e H(X, $) with ||*|| = 1 we have by (10.19) and (10.18) 

(i/2(^)-p'"/)^ > Im((i^2(^) -p'"/)^,*) 

> Im(|£|^-) -\y\snp\\W'{r + it)\\>\y\p 

te[0,y] 



2w-l 



where we have used that for y = the quadratic form of W{z) is real- valued. So the 
kernel of H2{r + iy) — p^"^ is trivial for y ^ 0. □ 



Lemma 10.6. For z e 7 and I e N 

(10.20) (^'^ 



j=0 



z-p 



j-i 
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where 



{2w)-\ 
1 ^ 

— T 



1 (-1 



E 



2w \ ( 2w 



giH Vqv=3 



2w 
qp + l 



i = 0; 

, i>o. 



The series in (10.20) converges absolutely. 
Proof. A striaghtforward calculation gives 

1 



2w „2w\~l 



- p 



;i0.21) 



) p2«,i ( Yl 



2w\ f z — p 



{2wy \ p 



q J \ p 

I 



If 2w e N, then the series on the right hand side is finite. Otherwise, by (10.17) and 
(10.6), for z e 7 the ratio of absolute values of any two sequential terms of the series 
satisfies 

z- p ( 2w \f 2w \ 
So, again by (10.17) and (10.6), we have 



z — p 


\2w-q-l\ 


p 


q + 2 



q>l. 



s—^ I 2w \ / z — p\i 



< 



2w 
2 



\Z — p\ X - 1 ^ Aw 



p 



q=0 



Thus we can decompose the expression on the right hand side of (10.21) into an absolutely 
converging series obtaining 



p=0 



l\ 1 



p J {2w)P V 5 \g + 1 



2w \ / z ~ p\i 



P 



{2wyi p^) + E(2^)p(p) ^^5]; (g, + i)---(g^ + i)) 



which finishes the proof. 



9l,--,9p>l 
91 H \-qp=j 



2w 
^qp + lj J' 



□ 



Let S{z) := H2{z) - z"^""! in H(X, *). Then by (10.9) on 7 the symbol of S{z) admits 
the representaion 

(10.22) s(^) = ^ r]z^^-"(2{ai,^)+z-'{\ai\^ + \y.\^)y + w-^{z). 

v=l ^ ' 
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Relations (10.22), (10.18), (7.1), Lemma 5.1(ii), and (2.20) imply that everywhere 
inside and on 7 

(10.23) i^'^^^^i ~ -'+""'+^-'+°+, / > 0. 

A version of the Jacobi's formula states that for any differentiable invertible matrix- 
valued function F{^z) we have 

tr {F\z)F-^{z)\ = (det [Fiz)]^ (^det [^(2)])"' 

(it can be proved, for example, using the expansion of the determinant along rows and 
the induction in the size of F). 

Then by Lemma 10.6 and the residue theorem 

(10.24) 

J2 r(p;X,$)^+^ 

= ^ j(^^+'(det [if2(-2) -p'"/])'(det [H2iz) - p^"" l]y' dz 
= ^ j( tr [z'^'-'Hl.iz) {H,{z) - p'-l) dz 

^ ^ (/ tr [(2«;^2"'+^7 + z^'+'S'iz)) V (-l)'>5'(z)(^^"' - p'''")-'-']dz 

^ _L ^ tr [(2wz2"'+^7 + z^+^>S'(^)) 

00 00 

Y^{z- pY'-' J](-1)'+M,+,+,,p^+'-^-(^+'+^-)5'+^(.)] dz 



X 

i=— 00 j=0 

E ^ ^ [(2^r-+^/ + r^+^5'(r)) E(-l)^+Mi+.+,,p^+'--(^+'+^)5^+^(r) 
z=o ■ j=0 



r=p 



We can restrict the summation on the RHS of (10.24) to 

1+ j ^h.= {M + K + d+l + {d-l)ad^i-2w)/{l- a^+i). 

Indeed, using the trivial fact that for any linear operator A in the finite dimensional 
Hilbert space spanned by with d e T^{h>) 

|trA|^ II A|| card T2j(i/), 

estimate (10.23), and relation (2.20) we can see that the sum of the terms in (10.24) with 
l+j > Iq contributes only to the order 0{p~'^~^'^^^'^) in (10.14), and thus after integration 
in $ the corresponding term can be included into the remainder Rj^ of Section 6. 
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Formula (10.12) shows that in order to compute the contribution to the density of 
states from S(QJ)p, we need to integrate the RHS of (10.24) against du and We are 
going to integrate against d# first: 



/ d$ / (cardT<^j(i/)) ' V r(p;X,$) 
= / du{c^rdT^{u))-'Y.Y. 



K+1 



1+l+jjP 



l+l-2w{l+l+j) 



X tr 



(10.25) 



dr'- 



1=0 j=0 



0{p, 



[ d*(2«;r'"'+^7 + r^+'>S'(r)),S'+^(r) 
du 



-M+2w~d 



) + 



-TT 

n(,j)CardY„(i.)^^ 



r=p 
1)1+ J 



0{P\ 



-M+2w~d\ 



l\ 



l+l-2w{l+l+j) 



X tr 



f2«;r2-+^/ 5'+^(r)d*-^^:t^ / 5'+^+i(r)d* 
+ ^itt(i / S^^^^Hr)d^ 



r=p 



We will prove that the integrand of the exterior integral in (10.25) is a convergent 
scries of products of powers of p and In p. The coefficients in front of all terms will be 
bounded functions of X, so afterwards we will just integrate these coefficients to obtain 
the desired asymptotic expansion. 

Let us discuss, how S{r) depends on p, X and In order to do this, we first look 
again at (9.12). As follows from Remark 10.2, the product eg'^cpg'^ does not depend on r 
and and by (8.4) 



10.26) 



|V^ee''^0''||L^(Md) <Pn'''^- 



For any rj G &s+i the application of the finite difference operator V^, to a polynomial 
decreases its degree by 1. Hence formula (10.5) ensures that 



;i0.27) 



K+l 



i^^j i=v rii,...,nx+i^O a=l 
niH \-nK+i^i 

jiH hjd'Sii-'^ 

Here C,-^V.'.'.-!f„ ^(X; 6) depend on the coefficients of (9.13) and satisfy a uniform estimate 



ini-'-Hif+i 



\^ini---nK+ 



(X;0)| < jd){o'm+i+0+) _ 
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Now 

(v^)(e) = E (v^.)(c + E%) (v-^(| • +e\i- - 1 ■ 

0=0 p=l 

The factors (V'^eeife) satisfy the estimate (10.26). For r] G @s+i we have 

(v,(i-+0ir-i-r)"')^o 

= (II + + 0ir - 1€ + ^ir)"'(ii + ^ir - \&''G{t e, r^i 

where 

^(C; 0, 77) := II + e\'^ - |^|^» - II + 77 + 0|^- + II + 771^"' 

= -2«;(r7, 0)111^--^ 

+ E l^lc""'' ((2(1' ^) + l^n' - (2(1, r7 + 0) + |77 + e\-y + (2(|, rj) + . 

j=2 ^ 

In analogy to (10.27) we have 

oo K+1 

(^^5(7^, 77)) (I) = E E C,...n.,,(X;0,77)r2-^-^ n(«i^ 

1=1^ ni,...,nif+i^O a=l 
niH 

with 

(10.28) |Q,...n.+, (X; 0, 77) I < 

Altogether, applying relations (10.26) - (10.28) to (9.11) and (9.12) we obtain 
(10.29) 

k—l s s—^i oo 

/ J / J / J / J / J / J npiio—isni—riK+iy i > 

LQ,...,Ls&J i^=^ 'ni^ - -Vs+tL^®s+i P=^ i=0 ni,...,n/f+i>0 

K+l s+n 

^ ^(2.-2);.+.o+-+^.-P-i -Q (sin $„)"" J] (|£ + 17, + 0,1^- - 1^ + 77,1^-) "\ 

a=l v=l 

where 

\^m-Vs+^ffi-0s+^ (v.n\\ < i(a^+i+0+)-(s-|U-p)/3+0+ 
\ siJ,pii-o'''i-sni---nK^i\ '> J\ rn 
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According to Lemma 10.1, 



(10.30) 



^r'-'-{2{te,) + 2{r)^,d,) + \e,\')-' 

oo K+l 

i=0 ni,...,nx+i^O a=l 

niH hriK+i^i 



and here 



If now subsitute (10.30) to (10.29), we obtain 
(10.31) 

k—l s s—n oo 

~ E E E EE E ^s/ipi(,o^-'^tsni-nif+i(-^' ^' 

a=l t)=l 

with 

|^'7l-^,+M^l-*'«+M fX-fl Tl ft W < J{am+l+0+)-(s-IJ.-p)l3+0+ 

The first sum in (10.22) can be written in the form 

$:(;;■' .-"(2(a,#) + .-'(|ar + |Xr))" 

D=l ^ ^ 

(10.32) oo K+i 

i=0 ni,...,nx+i>0 a=l 
niH hTiK+i<«+l 



where 
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Substituting (10.31) and (10.32) into (10.22) we can calculate the series for the symbol 
of the operator for / e N: 

(10.33) 

0u-,Bfe&^ 9=1 

/ j/fc h—u h—v—fi 00 

= EE E E E E E E ^(x;''....) 

iy=0 /i=i.^^_ ^^gjM=0 </>!,. ..,0ft_^+^e©^ p=0 i=0 ni,...,nK+i>0 

0i,...,6lh_^+^e©^ ni+-+nK+i<2^+p+j+/-i/ 

K+l h-u+fj. 

X ^(2-2^)^+(2^-l)/-.+.i+...+.,-p-. -Q (gin^^)n„ J-J ^2(e, 6>„) + 2((/>„, 6>„) + |6I,|2)"\ 

a=l v=l 

with 

|C(X-0 )| < p(/-''+»)("m+i+0+)-(/i-i/-;i-p)/3+0+ 

Note that the last product on the right hand side of (10.33) is of the form 



Here we have expanded the inner products d^) using Lemma 7.5(ii). The coefficients 
{6*} in the decomposition {dy)<ijx — Yliq^\i^q same sign and satisfy (7.4). 

Without loss of generality we may assume that all 6* are non-negative. The numbers 

satisfy p^^+Vn" Pn'^+'pV', since 

This numbers depend on X, but not on $ or p. The numbers kt — k{b\, . . . , 
positive, integer, and independent of ^. 

The following lemma is identical to Lemma 10.4 of [5], where for our purposes we have 
replaced the explicit constants 1/2 and 2/3 by -& and respectively. 

Lemma 10.7. For 1 ^ K d - 1; rii, . . . , uk+i e Nq; ki, . . . , kr e let Q := Ym=i ^t, 

(sin$i)"i . . . (sin<l)i^)"^(sin<l>A-+i)"^+i 



J, 



K 



-I 
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Then there exist positive numbers 60, px, and qx depending only on the constants (2.20) 
and K such that 

K 00 

where 

These estimates are uniform in the following regions of variables: 

(/ <lf< / o~^° < 6* < 0*^° o^"''^ < o 

Now using Lemma 10.7 we can compute the integrals of (10.33) over the domain 
G Mp} (recall that this integration is not needed for K = Q hy Remark 10.4). 
Substituting the result into (10.25), integrating in dv over f2(5J), and taking into account 
(10.12) and (10.14) we obtain in the region 2p„/3 < p < 6pn 

vo\A;{p)-yo\A;{p) 

K {«o+l)fc 00 

-EE E E^^^/^7v™^'-'"'^'-'"-'----'"-^(inp)«+o(p,-^+2-'^), 

9=0 h=0 ,^gjj=0 

with the coefficients satisfying 

I'-'qhj I ~ Pn 

This, together with equations (10.3), (6.9), Lemma 6.1, relation (2.20), Section 11 of 
[5], and the observation that the number of different quasi-lattice subspaces 2J is < pj'j^, 
completes the proof of Lemma 3.1 and, thus, of our main theorem in the case oi B = B 
with the symbol satisfying (2.8). As explained at the end of Section 3, the summation 
over J may be replaced by summation over Jq. 

It remains to relax the assumptions on B. This will be done in the subsequent section. 

11. Approximation 

In this section we prove Lemma 3.1 and thus Theorem 2.6 for general B using the fact 
that the proof is complete for B whose symbol fulfills the extra assumption (2.8). 

1. Given B satisfying the hypothesis of Theorem 2.6 and the number M, we fix the 
values of k and k in such a way that Lemma 3.1 holds true ior H = (— A)"' + S, where 
the symbol b oi B satisfying (2.8) is constructed at the end of Section 3. For > let 
us define (recall (6.1)) 
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We start by estimating the quadratic form oi B — B. For any G H^"'(]R'^) 

By Condition (2.14), the symbol of {B - B)7%^ satisfies 
Now Propositions 4.1 and 4.2 imply that 

(11.2) \\{-^ + l)--l\B-B){-^ + l)--l'7'^^\\^Cp-K 
Hence 

\{i:,9u,{B-B)T^^l^)\ 

(11.3) = |((-A + 1)-/V, J'«o(-A + 1)--/^(S - 5)y^„(-A + l)-/^(-A + 1)-/V) I 

^Cp-'=(V',(-A + 1)-/V), 
and the analogous estimates hold for the last two terms in (11.1). Thus (11.1) implies 

(11.4) \B-B\^B^''\ 
where B^'^'> is the operator of multiplication by the function 

(11.5) «'-\cp-'=(l+la'"^ i«i>«. 

in the momentum space. 

In view of Lemma 2.3(a), wc conclude that 

(11.6) Ar((-A)"' + 5, A) ^ Ar((-A)"' + S ± B^^\ A). 

So to prove (3.1) it will be sufficient to show that for p & In (which we assume everywhere 
below) the right hand side of (11.6) does not differ from A^((— A)"' + B,p^'^^ by more 
than 0{p-^). By (3.1) and Remark 2.7, it is enough to prove that 

(11.7) Ar((-A)"' + B± B^^\ A) = N{{-/^Y + B,X + Oip"^"^-^)). 
2. We note that for 

(11.8) := (4pf ^)V(— ) 
we have 

(11.9) N{{-AY + B±B^''\\) = iV((-A)"' + 5±yi^„sW±y^^5(^), A + 0(p2«'-")). 
Indeed, 

in view of (3.10). 
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3. Now we are going to prove that 

A^((-A)^ + B± ^ijoSW ± A) = A^((-A)^ + B± J^^^S^, A + 0{p''^-'^-^)). 

This will be done with the help of the following lemma, which is a development of 
Lemma 3.1 from [3]. 

Lemma 11.1. Let Hq, V , A he pseudo-differential operators with almost-periodic coef- 
ficients. Suppose that H := Hq + V is elliptic, selfadjoint and bounded below, and there 
exists a collection of orthogonal projections {-P/}^o commuting with Hq such that 

L 

(11.11) J2Pi = I arid Vni := PnVPi = for \l - n\ > 1. 

1=0 

Suppose that A — PqA and that 



a := \\A\\ < oo. 

At last, suppose that for A e R 

(11.12) Di -.^ dist {X, a{PiH Pi)) -{4 + 2^-^)a>0, Z = 0, . . . , L - 1 

and 

(11.13) max {a+\\Vii_i\\ + \\Vii+i\\)/Di^l/4. 
Then for 

(11.14) £:=2^-'^a 

we have 

(11.15) N{H, X-e)^ N{H + A,X)^ N{H, X + e). 

Proof. We will prove the first inequality; the second follows by interchanging the roles 
of Hq and + A. Let E\ be the spectral projection of (— oo, A] for H. By Lemma 4.1 
of [5] it is enough to prove that 

(11.16) (0,(i/ + A)0) ^ A||0||2 for every cj) e Ex-M^")- 
Let 

(11.17) S := min{a, 2'^'^ min A}, K := [2a/S] + 2, 

Q^I^Zj — 1 

so that 

(11.18) 2a^{K-l)d^3a 
and by (11.13) 

(11.19) K -1^ 3a/5 ^ 3max{l,2^+^a/ min A} ^ 2^+^ 

0^1<L-1 
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For (p G E\-^ introduce 

4>^ := {Ex-e-{k-l)S - Ex-e-ks)(p, k = l,...,K -1, 

K-l 

0^ := Ex-e-iK-i)s<f>, 0' := - 0^ = 5^ 0^ 

fe=l 

Then = X)f=i ®0'' and, letting 

(11.20) 77^= (A -£- (A; -1)5)/, k^l,...,K -1, 
we have 

(11.21) l|rylK5||</>1|. 

Let P-i :— Pl+1 ■= 0. Projecting (11.20) with Pi we obtain 

= Vn-i(f>U + {PiHPi -{\-s-{k- l)5))0f + Vii+^cPl,, / = 0, . . . , L, 
and thus by (11.13), (11.17) and (11.21) 

m\\ ^ iuw + ii^.-iiiii0tiii + \\vuMii\\)/Di 

^ 2-'-^U'\\ + J/4 + ||0t,||/4, Z = 0, . . . , L - 1. 

By induction, starting from Z = we obtain 

II ^ 2-2-^110^11 + 3||0t J/8, Z = 0, . . . , L - 1. 

Again by induction, using that ||0^|| ^ 110*^11, we get ||0f || ^ 2'--^||0'=||, Z = 1, . . . , L and 
thus ||0g|| ^ 2-^110^=11. Therefore, for A; = 1, . . . , X - 1, 

11^0^=11 = ||A0o^||^ 2-^0110^=11, 

and thus 

K-l 

P0'|| ^ ^ II 2-^v/F^a||0'|| 



k=l 

and 



K-l 

l'\\2 



|(0',^0')| = I Yl {<f>o,A^o) < 2-2^(X- l)a| 

k,m=l 

Hence 

(0, {H + ^)0) = (0', i70') + (0', A0') + 25R(0^, A0') + (0^, //0^) + (0^, A0^) 



^ (A - £)||0'||2 + 2-'^{K - l)a||0'|p + 2^-^ V:^^a||0'|| 110^11 
+ (A-£-(ii:-l)(^)||0^||2 + a||0^|p 

^ (A - £ + 2i-2^(X - l)a) ||0'||2 + (A - £ - (X - 1)5 + 2a) ||0-^||2 
^A||0f, 

where the last inequality follows from (11.18) and (11.19). □ 
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We now want to apply Lemma 11.1 to 

-.^ {-A)"" ±Tji^B^^\ V:=B, := ±7r^B'^''\ 

Note that 

(11-22) a := ||6||L^(Mdxs«,) + Mt^iwd,,^^^) 

does not depend on p„. For 

(11.23) L := [4 + log2 a + {M + d-2w) logs Pn] + 1 
we let 

(11.24) Rr.^Ro + lp'J", l^O,...,L-l, 
and introduce a family of projections 

(11.25) Po:^7r„ Pi:^'?R,--?R,_„ Z = 1,...,L-1, -.^ Tr^_^. 

Let us check that the hypothesis of Lemma 11.1 is satisfied. Relation (11.11) follows 
from (2.15) and (11.24). It follows from (2.9) that for Z ^ L - 1 

(11.26) WPiHPiW < \\PL_,{{-Ar + B)Pl-i\\ ^ 2||Pi_i(-A)'"PL-i|| ^ 2{RL_,f^. 
Also, for / ^ L - 1 

(11-27) + ||V;,+i|| ^ 2{Rl-i + pTf"". 

Since by (11.23) and (11.24) we have 

(11.28) Rl-i = i?o + [4 + log2 a + (M + d-2w) log^ Pn\pT ^ pT ^ogPn, 

relations (11.12) and (11.13) follow from (11.26) and (11.27) if pn is big enough. 
Applying Lemma 11.1, we get (11.15) with 

which implies (11.10). 

4. It remains to prove that 

(11.29) Nd-A)"" + B±7%^B^''\X) = N{{-A)'" + B,X + 0{pl'"-'^-^^)). 
Choose 

(11.30) s := p;'^-^. 
In view of (2.9), we have 

(-A)- + B^7r^{1± e){{-Ar + BpR^ © ^^^(1 ± l/e){{-Ar + B)Tr^. 
Therefore, 

3 (-A)- + B± S« ^ 7r. ((-A)- ± e{-Af + (1 ± e)BpR^ 

® ((-A)- ± {-Af/e + (1 ± l/e)B ± 
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Using (2.9) again and recalling the definitions (11.30), (11-5) and (11.8), we can estimate 
the last term on the right hand side of (11.31) from below: 

7%^ ((-A)"' ± {-Af/e + (1 ± l/e)B ± 5^)7^^ 

> ((-A)- - 2i-A)yep%^ ^ {Rl--2Rl^/e)7%^ ^ (Spn)'"^^., 

so it does not contribute to the density of states for p E In- For the first term we have 

((-A)- ± e{-Af + (1 ± e)BpR^ ^ 7r^{1 ± £)((-A)- + B)7r,, 

so 

A^((-A)^ + 5± y^^S^^^^A) 
(11-32) ^ N(pn^ ((-A)- ± ei-Af + (1 ± e)B) |,,^^l.(r<^)' ^) 

^^(3'«.((-A)'^ + ^)la>,.L.(M.),A/(l±e)), 

and the same estimates hold true for S'^'^^ replaced by 0. Combining these two versions 
of (11.32), we obtain 

(11.33) 

A^((-A)^ + B,X)^ Ar((-A)^ + B^ 'y'k^B^^\ X) 

^ iv(lP«. ((-A)- + B) |3,^^L.(R^)' V(l T e)) ^ N{{-Ar + 5, (1 ± e)\/{l t e)). 

Recalling that \ = p^"" ^ (4pn)2"' and (11.30), we arrive at (11.29). 
Combining (11.9), (11.10) and (11.29), we get (11.7). 
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